Electrophoresis of a polyelectrolyte through a nanopore 
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A hydrodynamic model for determining the electrophoretic speed of a polyelectrolyte through 
a nanopore is presented. It is assumed that the speed is determined by a balance of electrical 
and viscous forces arising from within the pore and that classical continuum electrostatics and 
hydrodynamics may be considered applicable. An explicit formula for the translocation speed as a 
function of the pore geometry and other physical parameters is obtained and is shown to be consistent 
with experimental measurements on DNA translocation through nanopores in silicon membranes. 
Experiments also show a weak dependence of the translocation speed on polymer length that is not 
accounted for by the present model. It is hypothesized that this is due to secondary effects that are 
neglected here. 



Electrically driven translocation of DNA across natu- 
ral and artificial nanopores can be detected on the single 
molecule level by observing the increase of electrical resis- 
tance of the pore during such events . Nanoscale pores 
may be fabricated by the self assembly of the natural 
protein a-hemolysin on a lipid bilayer membrane 0, Q 
or by annealing a microfabricated hole in a Si/SiOj sub- 
strate using an intense electron beam from a TEM 
The set up for detecting the translocation events consists 
of a reservoir containing an electrolyte that is partitioned 
into two chambers by a membrane with the nanometer 
scale pore forming the only communication between the 
two sides. An electrical potential difference is applied 
across the membrane and the resulting current is moni- 
tored. The passage of a DNA strand is signaled by a drop 
in the current. The duration and amount of these dips 
in the current contain signatures of the DNA such as its 
length and base sequence. Possible applications of the 
technique to rapid DNA sequencing is being explored Q . 

Lubensky and Nelson provided an interpretation 
of some of the features of the experimental work cited 
above. In particular, they proposed a drift diffusion equa- 
tion for P{s,t): the probability that the DNA is found 
to have a length s on a given side of the partition at 
time t. This equation was shown to explain qualitatively 
the shape of the observed distribution of translocation 
times. It was presumed that the drift velocity itself is 
determined by a resistive force that opposes the electri- 
cal driving force, but the exact physical nature of this re- 
sistive force remained unclear. Viscous resistance could 
be a reasonable contender for a resistive force localized 
around the nanopore, but a simple estimate appeared to 
indicate that its value was orders of magnitude smaller 
than what was required. 

Polymer translocation across nanopores driven by a 
variety of physical mechanisms have a number of other 
applications in biology; the injection of DNA from a virus 
into a host cell is a particularly interesting example 
These problems have been addressed by a number of au- 
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FIG. 1: Translocation of a polyelectrolyte across a nanopore; 
geometry of the pore region. 



thors lllllillllilllllllby formulating the prob- 
lem in a probabilistic setting as a transition between two 
states separated by a barrier in the configurational en- 
tropy of the polymer chain. In all cases, the hydrody- 
namic resistance of the pore when considered at all, is 
simply parametrized by a resistance coefficient. Most of 
the theoretical investigations on the subject to date are 
devoted to understanding how the translocation time of a 
polymer scales with polymer length, the scaling exponent 
being a quantity that may be expected to be indepen- 
dent of details such as the pore resistance and therefore 
amenable to experimental verification. 

In this paper, we explicitly determine the translocation 
time by actually solving a simple hydrodynamic model 
for the translocation of a polymer through a water filled 
pore. Translocation speeds are calculated and compared 
with the experimental data of Storm et al. for solid 
state nanopores. 

Figure ^ is a sketch of the set up that also helps to 
explain our notation. The pore is assumed to be cylin- 
drically symmetric j30| | about the z-axis and described 
by a function R{z) that gives the distance to the wall of 
the pore from the axis at any z-location. The part of the 
polyelectrolyte within the pore is modeled as a straight 
cylinder of radius a and carrying a uniform charge of den- 
sity A per unit length. The length of the pore is L, and 
the electric field, E{z), is in the z-direction. Constancy 



2 



of the current / through any cross-section orthogonal to 
the z-axis requires: 



E{z) = 



I 



aTr[R^z)-a^y 



(1) 



where a is the ionic conductivity of the buffer. It is as- 
sumed here that the electric field and current vectors 
make only a small angle to the z-axis, an assumption that 
is asymptotically approached in the limit |i?'(z)| ^ 1 for 
all z where the prime indicates derivative. The electrical 
conductivity a is taken as a constant in the pore with the 
same value as in the buffer. 

Next the fluid flow in the region between the poly- 
electrolyte cylinder and the walls of the nanopore is 
described, assuming that classical continuum constant 
density hydrodynamics is applicable. The assumption 
|i?'(z)| <C 1 immediately suggests that the problem may 
be treated in the lubrication limit [Ullil, so that (i) the 
pressure (p) is independent of the radial co-ordinate (r): 
p — p{z) (ii) the flow (u) is primarily along the z-axis, 
u ^ w{r, z)k and (iii) axial gradients are negligible in 
comparison with radial ones. Therefore, w satisfies 
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where p, is the dynamic viscosity of the buffer. Equation 
(0) does not contain an electric body force term because, 
except for a thin Debye layer next to the polyelectrolyte, 
the bulk solution is electrically neutral. 

In order to solve Q boundary conditions must now 
be specified. On the pore walls the classical boundary 
condition of no slip, is assumed: 



w{R{z),z) = 0. 



(3) 



The polyelectrolyte cylinder itself will be surrounded by a 
charge cloud of counter-ions. The thickness of this Debye 
layer is measured by the Debye length Ad, which for the 
high salt (~ 1 M KCl) buffers in these experiments is 
extremely small: Xjj ^ 0.3 nm. The problem is therefore 
treated in the thin Debye layer limit, which amounts to 
prescribing an apparent slip 0, 0| on charged surfaces 
(in SI units): 



w{a, z) 



<E{z) 



(4) 



Here v is the translocation velocity of the polyelectrolyte, 
C is the C-potential at the surface of the polyelectrolyte 
and e is the permittivity of the electrolyte. The boundary 
condition ^ provides the coupling between the fluid and 
the electrical problems. 

The solution to (O with boundary conditions (jSl and 
(El) is 
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where 



Tracer /i 



(6) 



defines a characteristic velocity scale for the problem. 
The pressure gradient p'{z) can be determined from the 
condition of mass conservation. 



w{r, z) (27rr) dr, 



(7) 



where Q is the flow rate through the pore. Substitution 
of © in lO gives 



dp 
Afi dz 



Q 2\nR^ 
^{Rl-l)f 



Ri - 1 



Rl - 21ni?* - 1 
{Rl - 1)/ ' 



(8) 



where f ^ Rl - I - {Rl + l)\nR^ and i?, = R{z)/a. 
The solution ^ and © still contains two undetermined 
parameters v and Q. These are determined by imposing 
the conditions that there is no pressure difference across 
the pore. 



p{0)^p{L) 



(9) 



and that the total force on the polyelectrolyte (which 
includes its Debye layer) is zero 







^ dw 
znap, -— 
or 



dz = 0. 



(10) 



The total force is zero since the polyelectrolyte moves 
through the pore without acceleration. 

Equations Q and ifTUIl yield after some algebra 
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(11) 

(12) 



where h, h, h-, h are dimensionless constants that de- 
pend solely on pore shape. They are defined as fol- 

lows: /o = (/-!>, h = {r\Rl - 1)>, h = {r\Ri - 

l)-\Rl - 2 In/?, - 1)) and h = {r\Rl - l)-\Rl - 
2 In/?, — 1)), where (•••)= L^^ /q (• • •) dz denotes aver- 
age along the pore length. 

TABLE I; Experimental parameters from Ref 



Ro 
(nm) 


a 

(nm) 


L 
(nm) 


ho 
(nm) 


a 
(deg) 


Ad 
(nm) 


A 

(e/nm) 


AV 

(mV) 


5.0/4.0 


1.0 


340 


40 


36 


0.3 


-5.9 


-120 



To calculate a numerical value for the translocation 
speed (u) from (|12|l , the ^-potential in © must be related 



to the charge per unit length of the polyelectrolyte, A 
which is usually known from the chemistry. In order to 
do so, one must adopt some model for the structure of 
the Debye layer. For the purpose of this calculation, the 
simplest model should suffice. Thus, it is assumed that 
the Debye layer is described by the Poisson Boltzmann 
equation in the Debye-Huckel limit Therefore, the 

electric potential due to the polyelectrolyte, 4>{r) satisfies 



ld_ 

r dr 



dr 



(13) 



with the boundary conditions 4>{a) — C, and 4>(oo) — 0. 
The solution to this boundary value problem may be ex- 
pressed in terms of the zero order modified Bessel func- 
tion Kq: 



(r) 



Koir/Xp) 
Koia/Xo) 



(14) 



The linear charge density A can be related to the potential 
(f) by Gauss's law: 



27ra ((>' (a) 



X 



(15) 



Evaluating (j)'{a) from (|14|l and substituting in (|15|l we 
have 



c = 



AAji Ko{a/XD) 
27rae ifi(a/AD) ' 



(16) 



where Ki is the modified Bessel function of order one. 
For the purpose of comparison with experimental data, 
it is convenient to replace the current / in © with the 
potential difference across the pore AV^ = ^(0) — V{L), 
where V{z) is the externally applied potential in the pore. 
Such a relation (Ohm's law) is readily obtained on in- 
tegrating between z — Q and z — L: Ijij^aa}) — 
IW/ihL) where h = {{Rl - 1)"^)- Therefore ® for u, 
may be written in an alternate form that does not involve 
the current or the conductivity. 
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The quantity Ueo = (e^//i)(Ay/i) has a very simple in- 
terpretation, it is the velocity with which a particle of 
any shape and surface potential C, will move if placed in 
an unbounded fluid medium and acted upon by the av- 
erage electric field that exists within the pore [23 • The 
funnel shape shown in Figure ^ may be assumed for the 
solid state nanopores: R{z) = Ra + (L ~ ha — z) tana if 
z < L — ho and R{z) = Rq otherwise, a being the semi- 
vertical angle of the cone. The integrals Iq,Ii,I2, I3 and 
I4 are then evaluated numerically. 

The prediction for the translocation speed v given by 
(|12|l is now compared with a set of measured values re- 
ported by Storm et al. 0, ■ Table HI summarizes the 
various parameters needed for such a comparison (e is the 




FIG. 2: Translocation velocity [v) as a function of the di- 
mensionless gap width {Ro/a — 1) in the case of 11.5 kbp 
dsDNA (squares), 48.5 kbp dsDNA (triangles) and 11.5 kbp 
circular dsDNA (circle). Lines are the predictions from the 
model (a) assuming counterion condensation in accordance 
with the Manning theory for a line charge (solid) (b) ignoring 
counterion condensation (dashed). 



charge on a proton). The pore radius Rq ~ 4.0 nm for the 
experiment with circular DNA (see Figure^! but 5.0 nm 
in all other cases. The dielectric constant e/eq — 80 and 
dynamic viscosity n = 8.91 x 10~^ Pa s for the electrolyte 
are taken as those of water. Equation (|16|l then gives [SJ 
C ~ — 56 mV. For highly charged polyelectrolytes such 
as D NA, it has been shown by Manning p3| and Oo- 
sawa ^2^ that some of the counterions condense on to 
the surface of the polyelectrolyte reducing its effective 
charge to Acft = X/qb- For DNA, at room temperature, 
in weak salt solutions (Ad 3> a) and away from bound- 
aries, the Manning factor qb ~ 4.2. However, imder the 
conditions of the experiment, X]j ^ a and further, the 
polyelectrolyte is in the vicinity of a dielectric/conductor 
interface which has an effect on qs H^. Thus, the value 
of ( is uncertain by perhaps as much as an order of mag- 
nitude. The characteristic velocity Ueo defined by 1)171) is 
UeO ~ 15.7 mm/s in the absence of counterion condensa- 
tion [qb = 1) but it is reduced to 3.7 mm/s if qb — 4.2 
is assumed. 

The lines in Figure [21 show the translocation velocity 
V calculated from (|12|l as a function of Ro/a ~ 1. The 
experimental data points are determined using Table I 
and Figure 11 of [l^. The translocation velocities are 
obtained by dividing the DNA length (11.5 kbp and 48.5 
kbp) by the measured mean translocation times. The un- 
certainty in the velocity corresponds to the spread in the 
translocation times. Four different modes can be distin- 
guished 14J for the DNA crossing, the two fundamental 
modes denoted as Type 1 and 2 and the mixed modes 
Type 12 and 21. In Type 1 the DNA passes through 
the nanopore without any folds. In the case of Type 2 
translocation of a linear chain or in the translocation of a 
circular chain, the DNA is folded in half so that its effec- 
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tive length is reduced to L/2 from L. In that case the two 
parallel DNA fibers that simultaneously thread the pore 
may be regarded as equivalent to a single cylindrical poly- 
electrolyte which has the same ^-potential but a larger 
radius, Ceff. The effective radius is determined by the 
requirement that the same fractional area of the pore is 
blocked. That is, the effective radius is acfr = \/2a ~ 1-41 
nm. Such a reduction has been used to calculate Ro/a—l 
for the Type 2 and circular DNA cases shown in Figure|2] 
An uncertainty of ±0.5 nm for the value of Rq is assumed 
due to errors such as the departure of the nanopore from 
a strictly circular shape and possible presence of hydra- 
tion shells on the surface. 

Lubensky and Nelson and Storm et al estimated 
the viscous force as fv ~ {2TTaLfiv)/{Ro — a) ~ 0.5 - 
5.0 pN, if u '--^ 1 — 10 mm/s. When compared with the 
electrical force fei ~ XAV ^ 113 pN, it appeared that 
viscous friction could not balance the electrical traction 
and therefore other mechanisms were needed to explain 
the observed translocation rates. Figure [21 demonstrates 
that this is not so. Two effects intervene to lower the 
estimate for the electrical force and raise the estimate 
for the viscous force (a) due to shielding by counter-ions 
the electrical force is reduced by as much as an order of 
magnitude (b) the change in the flow velocity takes 
place primarily across the thickness of a Debye layer, A^i 
and not over Rq — a, so that the actual viscous force is 
larger by a factor of approximately {Rq — a)/A_D ^ 10. 
In fact, in the thin Debye layer limit, the solution for 
the fluid flow within the Debye layer guarantees an exact 
cancellation of the electrical force with the viscous force 
on every surface element dS of the polyelectrolyte. Fig- 
ureElalso shows that the translocation velocity does have 
a weak dependence on polymer length. In a certain range 
of polymer l eng ths, the data can be flt by a power law: 
V ^ L^^-^^ |2J|. This fact cannot be explained by the 
model presented here or by any other model that localizes 
all of the resistive force within the pore. The dependence 
on polymer length must arise from additional effects not 
considered here such as the viscous resistance from the 
part of the polymer that lies outside the pore 2^ . 

The largest source of uncertainty in the above calcu- 



lation arises from the difficulty of obtaining accurately a 
value of C- The current understanding of the physics 
of the Debye layer is still incomplete so that even in 
the classical problems of electrokinetics, such as in the 
electrophoresis of a sphere, unresolved discrepancies ex- 
ist between Q and the C determined from more direct 
measurements of charge |2^. The use of the lubrica- 
tion equations and the assumption of a simplified axi- 
symmetric geometry for the pore also contribute to the 
error but these are likely to be much less than the ones 
just mentioned. One may question the use of continuum 
hydrodynamics in the first place to calculate the mean 
motion of the polymer. Note, however, in the case of 
water, the intermolecular spacing is of the order of 0.1 
nm which is significantly smaller than the ~ 2 nm inner 
diameter and 4 — 5 nm outer diameter of the pore. The 
applicability of the no slip condition at the solid liquid 
interfaces is still a matter of contention Neverthe- 
less, possible slip lengths are miniscule and amounts to 
an uncertainty in the values of Rq and a by perhaps a 
fraction of a nanometer. In fact, classical hydrodynamics 
(i.e. Stokes equations with no slip boundary conditions) 
work reasonably well for water down to several tenths 
of nanometers; for example, for non-electrolytes, molec- 
ular sizes calculated on the basis of the Stokes-Einstein 
relation or Einstein's viscosity law for dilute suspensions 
agree with molecular structure based determinations to 
within 10— 15 percent even for molecules in the 0.3- 
0.5 nm range. The use of the continuum hydrodynamic 
model for calculating statistical averages in the manner 
used here is therefore not likely to be a significant source 
of uncertainty, at least for the solid state nanopores. 

In summary, a model of the pore resistance based on 
continuum hydrodynamics and electrostatics produces 
estimates for the translocation speed of dsDNA in solid 
state nanopores to within an order of magnitude of exper- 
imental values. The present analysis needs to be mod- 
ified for protein nanopores because certain approxima- 
tions such as the thin Debye layer are not applicable in 
that case. 
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